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Wright 分 数 阶 时 灌 微 分 方程 的 离散 化 过 程 
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摘 要 : 引进 了 一 种 离散 化 方法 对 分 数 阶 时 滞 微 分 方程 进行 离散 化 求解 。 首 先 考察 Wright 分 数 阶 时 滞 微 分 方程 ; 其 次 
分 析 相应 具有 分 段 常数 变 元 的 Wright 分 数 阶 时 灌 微 分 方程 ， 并 应 用 离散 化 过 程 对 模型 进行 数值 求解 ; 然后 根据 不 动 点 
理论 讨论 该 合成 动力 系统 不 动 点 的 稳定 性 ; 最 后 借助 MATLAB 对 模型 进行 数值 仿真 ， 并 结合 Lyapunov 指数 、 相 图 、 
时 间 序 列 图 、 分 岔 图 探讨 模型 更 多 复杂 的 动力 学 现象 。 结 果 显 示 ， 提 出 方法 成 功 对 Wright 分 数 阶 时 灌 微 分 方程 进行 离 
pe 
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On discretization process of fractional-order delay Wright equation 


Liu Shanshan, Gao Fei 
(School of Sciences, Wuhan University of Technology, Wuhan 430070) 


Abstract: This paper introduced a discretization process to discretize fractional-order delay differential equations. First of all, 
it considered the fractional-order delay Wright fractional differential equation. Then, it analyzed the corresponding fractional- 
order Wright differential equation with piecewise constant arguments and applied the proposed discretization process to solve 
the model numerically. After that, according to the fixed points theory, this paper investigated the stability of the fixed points of 
the resultant dynamical system. Finally, it carried out a numerical simulation including Lyapunov exponent, phase diagrams, 
time series diagram, birurcation using Matlab to reveal more complex dynamics of the model. Simulation experiments show that 
this paper succeeds in discretizing fractional-order delay Wright equation. 
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0 az 节 、 农 业 商品 市 场 、 非 线性 光学 和 血液 中 的 噬 中 性 粒 细胞 。 我 
们 将 特别 关注 这 样 的 系统 ,其 中 时 混 微 分 方程 表达 为 如 下 形式 : 

混沌 被 认为 是 继 量子 力学 与 相对 论 之 后 的 第 3 大 科学 的 发 x (n) f(x().x(r-2)). x(t) eR",t>0 a) 

现 。 混 沌 系统 是 非 线性 动力 学 映射 主要 表现 形式 之 一 ， 具 有 良 这 里 在 时 间 t 时 系统 的 状态 为 x(1) ， 在 过 去 的 时 间 状 态 


好 的 类 随机 、 非 周期 、 对 初始 值 敏感 、 历 经 各 态 并 可 确定 等 特 KAH (t-r), RET 为 一 个 固定 的 时 间 延 迟 。 此 外 时 滞 微 
性 ， 被 信息 安全 、 扩 频 通 信 及 工业 控制 等 领域 的 研究 者 们 广泛 。 ”分 方程 出 现在 很 多 数学 模型 中 : 例如 种 群 动态 学 〔 妊 娠 时 间 
K ANTARE iam. JL HER, AMD CJESEX dg 。 ”考虑 在 内 )， 传 染病 (传染 周期 考虑 在 内 )， 生 理 和 药理 动力 


沌 系统 迅速 发 展 并 且 引 起 学 者 们 的 极 大 关注 。 分 数 阶 微 积分 被 ”学 〈 血 液 循环 周期 )， 化 学 动力 学 等 RHRSRRPUSR 。 

看 做 是 整数 阶 到 非 整数 阶 微 积 分 的 推广 ， 已 知 当 系 统 阶 数 超过 分 数 阶 微 积分 是 传统 意义 的 整数 阶 微 积分 的 延伸 和 拓展 ， 
3， 系 统 出 现 混沌 行为 。 但 研究 人 员 发 现 , 对 于 一 些 分 数 阶 微 积 。 是 将 求 导 次 数 和 求 积 次 数 由 整数 推广 到 非 整 数 。 这 里 “分 数 ” 是 
分 总 表现 出 阶 数 小 于 3 的 混沌 行为 , 如 Hartley S A ARRAS 一 个 广义 的 概念 ， 不 仅 包 含 分 数 ， 还 包含 无 理 数 和 复数 等 ， 所 
的 研究 前 述 了 该 现象 。 以 "分数 阶 微 积分 其实 是 “ 非 整数 微 积分 "。 分 数 阶 微分 是 分 数 


其 中 ， 时 滞 微 分 方程 DDES) 是 依赖 于 当前 和 过 去 历史 状 ” 阶 微 积分 重要 研究 领域 之 一 ， 也 是 处 理 实际 应 用 中 非 线性 问题 
态 的 动力 系统 的 数学 表达 式 ， 例 如 延 时 反馈 神经 系统 、 呼 吸 1 的 有 力 工具 ， 本 文 就 是 基于 分 数 阶 微分 的 基本 理论 展开 对 
Wright 时 滞 微 分 方程 非 线性 现象 的 研究 和 分 析 。 由 分 数 阶 微分 
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的 定义 可 以 看 出 ,分 数 阶 微分 
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与 以 前 所 有 采样 点 的 采样 值 有 关 ， 


取 值 有 关外 天 | 此 分 


记忆 、 遗 传 等 特性 的 过 程 ， 区 别 于 整数 阶 微分 系统 ， 也 是 研 


而 整数 阶 微分 只 与 当前 采样 点 和 当前 时 刻 的 前 一 时 刻 采样 点 


的 


前 ， 对 于 分 数 阶 微分 方程 的 求解 主要 是 采用 频 域 法 weeamsm 


9 和 时 域 法 * ssew 两 种 方法 。 但 近来 研究 发 现 相 比 第 一 种 


数 阶 微分 系统 适合 刻画 描述 具 


分 数 阶 微分 系统 的 必要 因素 。 


再 有 整数 阶 系统 是 分 数 阶 系统 


次 取 整 数 时 的 特例 ， 是 现实 系统 理想 化 的 一 种 处 理 方式 ， 所 
所 有 混沌 系统 都 可 以 说 是 分 数 阶 混沌 系统 。 因 此 采用 分 数 阶 


沌 系统 对 混沌 现象 进行 描述 更 具有 普 适 性 。 而 本 文 将 Wright 
消 微 分 方程 推广 到 分 数 阶 领域 ， 其 离散 化 过 程 对 分 数 阶 求解 
供 了 理论 指导 ， 为 分 数 阶 微 积 分 的 应 用 开拓 了 新 领域 。 近 几 年 
因为 分 数 阶 微 积分 在 科学 技术 领域 的 广泛 应 用 ， 吸 引 到 很 多 学 


有 
y 
阶 
以 
混 
时 


提 


Br 


学 


者 , 他 们 利用 数值 仿真 研究 该 类 方程 。 因为 Caputo 型 分 数 阶 微 


分 描述 问题 的 初始 条 件 与 整数 阶 微分 方程 是 一 致 的 ， 因 此 本 文 
] Caputo 型 分 数 阶 微分 进行 研究 ， 首 先 立 述 Caputo 型 分 数 


BE 


阶 微 积分 的 定义 : 


定义 1 eign 函数 (1) 定义 在 区 间 (0,+%) ， 其 中 


数 Be R 的 分 数 阶 积分 公式 为 


分 数 阶 微分 的 公式 为 


D'f(r)-1"D'f(t), D-— 


Xp e(n-Ln), nod 
文献 错误 ! 未 找到 引用 源 。 
x( 


d 
dt 
FE 整数 ，1 >0。 
研究 了 整数 阶 Wright 方程 


1)=-px(1—1)(1+x(7)) 


方法 ， 第 二 种 更 加 有 效 ， 因 第 一 种 方法 在 检测 
的 。 


二 < 


昆 沌 不 总 是 可 靠 


通常 分 数 阶 微分 方程 的 解析 解 不 总 是 合意 的 ， 于 是 不 得 不 
借助 数值 解 的 方法 并 辅助 计算 机 仿真 。 分 数 阶 微分 方程 数值 解 
的 方法 主要 包括 : 有 限 元 方法 、 变 分 迭代 法 、Adomain 分 解法 、 
预 估 -校正 法 等 。 其 中 在 文献 错误 ! 未 找到 引用 源 。 分 别 对 分 数 阶 
logistic 模型 和 分 数 阶 Henon 模型 给 出 了 分 数 阶 的 数值 求解 方 
法 , 但 该 方法 与 预 估 - 校 正法 类 似 , 求解 结果 相对 其 它 算法 相 比 
较为 精确 ， 但 是 两 者 都 存在 计算 结果 复杂 ， 并 且 前 者 的 数值 解 


不 仅 与 当前 状态 相关 而 且 依赖 过 去 所 有 状态 ， 故 在 进行 数值 仿 


阶 


Q) 


其 中 : p 为 实 参数 ， p» 0， 表明 方程 2) 在 w =F RE Hopf 分 


贫 ， 并 且 在 a > 二 时 存在 周期 解 。 文 献 错误 ! 未 找到 引用 源 。 证 


明了 Wright 微分 方程 可 由 Logistic 方程 经 过 一 系列 变换 而 来 ， 


但 Wright 微分 方程 与 Logistic 


模型 一 样 看 似 简单 


但 是 参数 P 在 一 定 范围 内 变化 时 ， 它 有 着 极为 复杂 的 动力 学 


为 ， 其 中 包括 了 分 贫 和 混沌 。 


本 文 基于 Caputo 型 分 数 阶 微分 定义 和 性 质 ，; 


且 可 确定 ， 


p 


fT 


各 整数 阶 


Wright 时 滞 微 分 方程 推广 到 分 数 阶 Wright 时 灌 微 分 方程 ， 此 
外 将 时 消 推 广 为 一 般 形 式 ， 改 后 的 方程 表达 式 如 下 : 


D'"x(t) e -ex(t - c)(1* x(t)).re (0,7] 


初始 条 件 为 


这 里 a e (0,1] 。 然 后 根据 动力 系统 不 动 点 稳定 性 定理 和 


混沌 动力 学 分 析 研 究 分 数 阶 Wright 时 灌 微 分 方程 的 动力 学 性 


质 以 及 混沌 现象 。 


1 ”离散 化 过 程 


G) 


x(0) 2 x, (4) 


当 分 数 阶 微分 模型 建立 2 


后 ， 便 面临 着 如 果 求 解 的 问题 。 


真 加 大 实验 难度 ， 进 而 较 难 分 析 模 型 的 动力 学 现象 。 


c 


MEFR, RA AAA BUB E 0A) 71 Fe 56 SUB 


越 多 生态 数学 学 者 的 关注 ， 这 一 模型 正解 的 稳定 性 、 有 界 性 、 


Wc PEL b EB] T tA D y Roo mios 
本 文 参考 文献 错误 ! 未 找到 引用 源 。~ 错 误 ! 未 找到 引用 源 。 


引用 源 。 


的 离 


散 化 过 程 , 定义 了 有 共有 分 段 常数 变 元 的 分 数 阶 Wright I itr) 


方程 ， 并 对 其 进行 离散 化 处 理 ， 方 程 表达 如 下 : 


D'x(r)- -ad 全 - Jos H 小 e(o.] 


初始 条 件 为 式 (4)。 
将 方程 (5) 逐 段 离散 的 过 程 如 下 : 


邻 1e[0,.7)， 那 么 <e[0,1) ， 故 有 

r 

D*x(t) 2 -gx, (-c)(1- x) 
那么 方程 (9) 的 解 为 


x(t) -2x-I* (ex (7c) (14 x,)) 


esiste ra 


t 


= Xo -a (os (Ot) 


T(1*« 


Ate[rar), 那么 -< [1.2) ， 故 有 
D*x(t)=-px (r-r)(1+x (r)) 
那么 方程 (6) 的 解 为 
x (i) =x (0) 12 (ox (c) (ox (7) 
-xi)-(ostr- sten (EET 


t-r) 


-x (r) -E {ox (r-l) 


(+a) 


重复 上 述 过 程 ， 易 得 到 


(6) 


(0) 


xa (t) =x, (nr)- l os (nr- r)(1 X, (nr))).t € [nr. (n +1) r) 


T(I+aw) 


4 t—(nl)r 
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X. ((n 十 1) r) (nr) 一 TET (ox, (nr 一 z)(1 tX, (nr))) 
4 7 三 六 ， 故 有 
Xa =X ls (+x, )) (8) 


值得 注意 的 是 欧 拉 离散 化 方法 是 对 求 导 的 近似 估计 ， 而 预 
估 - 校 正法 是 对 积分 的 近似 估计 和 中 3,'， 而 本 文 提出 的 离散 
化 方法 是 对 当前 研究 系统 区 间 右 端的 一 种 近似 逼近 。 此 外 ， 可 
以 注意 到 w 一 1 时 ， 即 欧 拉 离散 过 程 。 


2 Wright 时 滞 微 分 系统 的 动力 学 分 析 


对 于 非 线性 系统 ， 可 采用 在 相 空 间 观察 其 轨道 的 方法 。 非 
线性 系统 随 着 时 间 的 演变 ， 将 趋 于 维 数 比 原来 相 空 间 低 的 极限 
集合 -吸引 子 。 通 常 简单 的 吸引 子 有 不 动 点 、 极 限 环 和 环 面 。 简 
单 吸引 子 又 受 控制 参数 的 影响 ， 随 着 控制 参数 的 变化 ， 简 单 吸 
引子 可 逐渐 发 展 为 奇怪 吸引 子 ， 此 时 系统 是 混沌 的 并用。 
这 种 当 控制 参数 变化 到 某 个 临界 值 时 导致 系统 的 动力 学 状态 发 
生变 化 的 现象 称 为 分 岔 。 下 面 从 系统 不 动 点 的 稳定 性 开始 ， 考 
察 系统 不 动 点 稳定 性 ， 分 岔 点 位 置 ， 分 岔 类 型 以 及 发 生 分 贫 后 
的 混沌 判别 。 

接 下 来 研究 方程 (8) 的 定点 ， 通 过 求解 下 式 ， 可 知 其 存在 两 
个 定点 ， 分 别 为 0 和 -1。 


2i (ex 十 x)) 


xX-—x-— 


为 研究 系统 (8) 定 点 的 稳定 性 ， 给 出 如 下 定理 : 
定理 | "AUSRSUSUBM. 假设 函数 f 为 定义 在 R 上 的 一 个 光滑 
映射 ， 为 为 函数 三 的 的 定点 。 


(x) «1 


， 那 么 如是 系统 的 稳定 点 ; 


(x)|»1 
首先 考虑 定点 nl 的 情况 ， 


， 那 么 为 是 系统 的 不 稳定 点 。 


可 知 此 时 


f'x)- 吕 >1， 故 系统 在 定点 -1 是 不 稳定 的 ; 而 
| 


TON AN 
T(1+a) 
对 于 定点 "0 是 稳定 的 ， 当 且 仅 当 


2T(1+aw) 
r* 


0<p< (9) 


此 外 ， 为 数值 仿真 的 方便 ， 可 将 方程 (8) 改 写 为 如 下 二 
分 方程 : 


阶 差 


a 


py, (1+%,)) 


Aya 二 加 一 


FU (10) 


Ynn T Xn 

那么 此 时 系统 存在 两 个 定点 : jx =(0,0), fx =(1-1) 。 为 研 
究 定点 的 局 部 稳定 性 ， 需 求 出 系统 在 每 个 定点 的 雅克 比 矩 阵 的 
IE. DAR IR 2] 7) AR Gt Re E VE sg Sp RURHUUP RT 知 当 
| 外 <bi=1,2， 系 统 渐进 稳定 。 可 知 系统 在 任意 点 Gc y) 雅克 比 
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矩阵 的 表达 式 为 
PS D 1 
4 r(1«a) ^" rüsaj^ » (11) 
1 0 


2.1 定点 1 的 稳定 性 分 析 及 分 贫 
将 第 一 个 定点 fox, =(0,0) 带 入 便 可 得 到 


1 "m a 
J= r(1+a) a2) 
1 0 


得 到 的 特征 方程 为 


FU T (13) 
求 得 的 特征 值 为 4a 7050.5 j1- E E 。 可 知 当 


1-4 


r“p r(1+a) _ 
0 LR +0 2 <1 ~ 
Fü- >0 即 0<p< px 时 ， 此 时 0<4，。<1， 不 动 


a) 


点 为 稳定 结 点 ; Ton 中 时 ， 和 fill A, 


tou, HIN 


Al-ul- 


I9) 


可 知 当 pp 时 ， 此 时 |4|=1 ， 于 是 乘 子 开始 穿 过 夫 洛 


开 单 位 圆 ， 发 生 Neimark-Sacker 7 2i, 
rT(l+a) 


a 


JL PUEH, DAA 


A(x»)-(60). p= 


对 为 对 于 不 动 点 而 言 ， 当 满足 |%|<4i=1,2 的 充分 必要 条 
件 为 Jury Ja] se 2 (tp ries. ~ 错误 ! 未 找到 引用 源 。 表达 如 下 : 


a) F=1+T+D>0. 


b) TC=1-T+D>0. 
c) H=1-D>0. 
这 里 ， 刀 分 别 表示 矩阵 的 迹 和 雅克 比 行列 式 的 值 。 
T: Tr(J)2 Ja Js, 
D:Det(J) 2 JaJa- JaJa 


a 


式 (12) 可 知 ， 此 时 7 21, D 


r 
- Pues ， 可 知 此 时 条 件 


ajb) 总 是 满足 的 ， 而 对 于 条 件 c) 不 是 总 满足 的 。 根 据 Jury 判定 


准则 ， 当 1-D=0 时 ， 即 1-f(0)-1 


r? 
Fco)? ' 那么 


p=+ 人 9 ， 此 时 系统 存在 单位 元 上 一 对 共 辆 复 根 ,由 图 1.3 


1-4 
rur a) 时 ， 系 统 失 稳 发 生 Neimark-Sacker 7j 2 o 


可 知 当 p> 


2.2 定点 2 的 稳定 性 分 析 
第 二 个 定点 fos =(-1,-1) 的 雅克 比 和 矩阵 表达 如 下 : 
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1 
J- Via) j (14) 
1 0 
得 到 的 特征 方程 为 
Mis vn 


求 得 的 特征 值 为 : A-0, 7 


UTE Ps 可 知 无 论 e 取 任 


FÉ, EA fo -(7L-1) 总 
2.3 混沌 判别 

混沌 系统 在 相 空间 中 是 由 奇怪 吸引 子 的 不 规则 轨 线 来 描述 
PARRE, 为 了 定量 刻画 奇怪 吸引 子 的 运动 状态 ， 引 入 了 
Lyapunov 指数 。 


定义 1 错误 ! 未 找到 引用 源 。 


是 不 稳定 的 。 


IT 


般 二 维 映射 S W Jacobi ERE, 


设 表示 


^ 


I= F) np ona) - [^9]... 
将 J 的 2 个 复 特征 根 取 模 后 ,依次 从 小 到 大 顺序 排列 后 为 
lal 2|4] 


那么 f 的 Lyaounov 指数 定义 为 


À, - lim-log Aj? 


(k 212) (16) 


该 定义 是 计算 一 般 二 维 映 射 的 Lyaounov 指数 的 理论 基础 。 
Lyaounov 指数 的 大 小 表明 相 空 间 中 相近 轨道 的 平均 收敛 或 发 
散 的 指数 率 。 可知 ， 非 负 的 Lyapunov 指数 的 方向 ,对 吸引 子 起 
支撑 作用 , 而 负 的 Lyapunov 指数 对 应 着 收缩 方向 。 这 两 种 因素 
对 抗 的 结果 就 是 伸缩 与 折 著 操作， 这 就 形成 奇怪 吸引 子 的 空间 
几何 形状 。 因 此 ， 对 于 二 维 映射 的 奇怪 吸引 子 而 言 ， 其 最 大 
Lyapunov 指数 为 正 ; 反之 亦 然 。 下 节 数 值 仿真 给 出 在 不 同 的 
r,&, p FEES) Lyaounov 指数 ， 此 判定 系统 发 生 混沌 与 否 。 值 


«-0.95,7-0.25 
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ChinaXiv& fERBT! 
X]4444, 3: Wright Z^ Ip p Pie 427 77 42:89 A CICER 


得 注意 的 是 ， 对 于 系统 式 (8) HJ Lyaounov 指数 表达 为 


Lya.exp = lim > log, [ 一 ri^ 十 2) 
x alit, EBRA xu ox, rex (1x) 的 
Lyaounov 指数 的 表达 式 。 


3 ”数值 仿真 分 析 


本 章 主要 借助 MATLAB 进行 数值 仿真 来 说 明理 论 结果 以 
及 揭示 系统 (10) 在 z=7 复杂 的 动力 学 现象 。 这 里 主要 采用 控制 
变量 法 ， 即 在 其 余 参 量 保持 正常 值 的 条 件 下 单独 探讨 一 个 参量 
对 系统 的 影响 ss。 本文 探 讨 了 系统 参数 OLA 1~3 )， 
miir (图 4) 以 及 阶 数 & 对 分 数 阶 微分 的 影响 。 此 外 ， 一 个 系 
统 发 生 混沌 与 否 ， 一 个 关键 的 参考 依据 就 是 最 大 李 雅 普 诺 夫 指 
数 。 故 最 大 李 雅 普 诺 夫 指数 是 对 混沌 最 好 的 观测 者 ， 当 系统 发 
生 混 沌 时 ， 最 大 李 雅 普 诺 夫 为 正 ， 反 之 为 负 。 

对 于 系统 (10), 图 1 中 分 别 取 时 滞 上 = 0.25， 分 数 阶 阶 次 C 
取 定 为 0.80 到 0.95, 此 时 分 数 阶 阶 次 具有 一 般 性 ,改变 参数 P ， 
观察 在 此 区 间 内 系统 随 着 参数 改变 的 动力 学 现象 。 此 外 初 值 均 
设 定 为 =0.1, yo=0.2， 步 长 设 定 为 0.001。 可 知 随 着 参数 Pp 
的 增 大 ， 系 统 由 稳定 态 发 生 Neimark-Sacker 分 偏 ， 进 入 混沌 ， 
此 时 从 图 中 可 看 出 在 分 岔 点 Lyapunov 指数 为 0, 进入 混沌 后 大 
于 0。 此 外 对 于 很 多 动力 系统 而 言 分 岔 及 混沌 现象 都 是 由 时 清 
引起 的 ， 故 这 里 蔡 换 时 灌 为 r+=0.5,r =1， 改 变 参 数 pP ， 查 看 时 
灌 对 系统 的 影响 ， 见 图 2 和 3。 从 图 1-3 分 别 可 以 看 出 ， 随 着 了 
和 阶 数 & 的 增 大 , 系统 (10) 趋 于 稳定 , 不 同 的 是 + 增 大 ,系统 稳 
定 区 域 变 小 ， 分 岔 提 前 ， 而 阶 数 wx 反之 。 

从 图 1~3 可 以 得 到 分 数 阶 阶 次 的 改变 对 于 影响 着 系统 (10) 
的 运动 状态 , 先 初 值 设 定 为 加 =0.1yo=0.2 , 步 长 设 定 为 0.005。 
固定 参数 P ， 之 后 分 别 取 时 兆 了 为 0.14~0.3 Ere 
0.60~0.95， 得 到 分 贫 图 如 图 4 所 示 ， 可 知 随 着 c 的 增加 系统 趋 
于 稳定 。 
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a=0.85,r=0.25 a=0.8,r=0.25 


图 1 r=0.25 时 ， 分 数 阶 & 取 不 同 值 时 系统 (10) r 


0=0.95,r=0.5 0-0.90,7-0.5 
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a=0.85,r=1 


r=0.14, p=3.8 


当 参 数 变 化 到 某 个 值 时 , 系统 相 图 的 拓扑 结构 也 随 之 变化 ， 

在 图 形 上 表现 为 参数 取 改 值 时 相 图 发 生 突 变 。 接 下 来 给 出 系统 

(10 Æp=4, r=0.25, a=0.95 时 的 相 图 。 从 图 1 可 知 当 

P=4, 系统 为 混沌 状态 ,Lyapunov 大 于 0, 故 图 5 表现 为 混沌 
吸引 子 。 


图 3 r=1 时， 分数 阶 oc 取 不 同 值 时 系统 (10) 的 分 贫 图 


44, F: Wright 分 数 阶 时 灌 微 分 方程 的 离 


co=0.8.r=1 


0.95 0.955 0.96 0965 0.97 
a 


图 4 了 固定 ， 分 数 阶 & 变量 时 系统 (10) 的 分 岔 图 


r=0.25, o=0.95,p=4 
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5 r=025，D=4，w=0.95 时 系统 (10) 的 相 图 


混沌 动力 系统 最 重要 的 特性 之 一 是 对 初始 条 件 变化 十 分 敏 
感 ， 这 意味 着 对 当前 轨迹 的 任意 小 的 扰动 都 可 能 被 放大 ， 并 在 
多 次 迭代 之 后 导致 显著 不 同 的 未 来 行为 。 从 


(olano (xy) FERRAR, MEEN 
; 


日 欧式 距离 衡量 ， 表 达 如 下 : 


3 2 
Ed, „ = (s -x) «(» -x) (17) 


为 了 显示 混沌 系统 对 初始 条 件 的 敏感 性 ， 图 6 Za 


两 个 无 


限 接近 但 不 相同 的 初始 点 (0(0.1,0.1),0"(0.101,0.1)) 之 间 的 距离 


(用 红线 表示 )。 


0.03 r —r 


RARA 
一 稳定 系统 


0.025 上 


0.02 凡人 


5 0015 上 


0.01 上 


0.005 
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"s 


6 系统 (10) 的 时 序 图 


显然 ， 对 于 混沌 系统 而 言 ， 尽 管 对 初 值 扰动 很 小 但 在 几 次 
迭代 后 被 放大 。 而 由 稳定 系统 生成 的 另 一 蓝 色 基 准 线 ， 尽 管 两 
条 轨迹 从 两 个 差异 略 大 的 点 o(0.1,0.1) 41 '(0.1,0.12) 出发 但 
随 着 迭代 进行 ， 两 条 轨迹 之 间 的 距离 趋 于 0， 表 现 出 明显 的 初 
台 敏感 性 ， 进 而 验证 了 上 述 的 结论 。 
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4 ”结束 语 


本 文 研究 了 分 数 阶 时 沾 Wright 方程 , 并 运用 了 一 种 离散 化 
方法 对 其 进行 离散 求解 。 并 以 & ， 为 参数 通过 数值 仿真 ， 并 
通过 分 岔 图 、Lyapunove 指数 、 相 图 对 其 动力 学 
另外 注意 到 ， 当 w -1 时 ， 即 欧 拉 离散 过 程 和 fs 。 但 是 区 

别 于 欧 拉 法 的 是 ， 本 文成 功 对 二 阶 Wright 微分 方程 进行 离散 。 
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